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ЗАДАЧА ТЕРРАУПРУГОСТИ В ПОЧВООБРАБОТКЕ 


Решается задача построения математической модели процесса взаимодействия ра- 
бочего органа с почвой с учётом упругости его крепления. Отмечается аналогия её 
задачам аэро- и гидроупругости в классической механике. Модель построена в виде 
векторно-матричных уравнений в форме Коши, для чего введено понятие вектора 
состояния упругозакреплённого рабочего органа. Приведены перспективы при- 
менения модели для описания динамических эффектов почвообработки. 
Ключевые слова: почвообработка, рабочий орган, упругость, динамика. 


Введение. Применяемые в почвообработке математические описания вза- 
имодействия рабочего органа с почвой широко используют постулат, что 
рабочий орган закреплён жёстко, движется поступательно, не отклоняясь 
от прямолинейной траектории [1]. Это расходится с реальным процессом. 
Исследования показывают, что движение рабочего органа в почве сопро- 
вождается значительными упругими смещениями и колебаниями интенсив- 
ностью до 89 по виброускорению [2]. Упругие смещения соизмеримы с глу- 
биной хода и оказывают сильное влияние на качественные и энергетиче- 
ские показатели почвообработки. Величина упругих смещений составляет 
от 10 до 110 мм; искажения геометрии рабочих органов достигают 
20 град/кНн [3]. Абсолютно жесткого крепления не существует, поэтому в 
строгом смысле все крепления рабочих органов почвообрабатывающих ма- 
шин упругие, и взаимодействие рабочих органов с почвой в земледельче- 
ской механике должно рассматриваться обязательно с учётом упругих сме- 
щений. 

Новый подход к моделированию процесса взаимодействия рабочего 
органа с почвой связан с учётом его упругих смещений. Это позволяет в 
почвообработке рассматривать взаимодействие рабочего органа с почвен- 
ной средой как задачу терраупругости по аналогии с задачами аэроупруго- 
сти и гидроупругости в механике, но для среды с более сложными реологи- 
ческими свойствами. 
Постановка задачи. Расчётная схема взаимодействия упругозакреплён- 
ного рабочего органа (УРЗО) с почвенной средой представлена на рис. 1. 
Рабочий орган 3 в виде трёхгранного клина крепится к раме 1 консольной 
стержневой системой 2 разной степени сложности (стойка, подвеска, гря- 
диль и т.п.), выносящей место крепления за почвенный горизонт. Рама 
движется поступательно со скоростью \. 
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Рис.1. Расчётная схема упругих смещений рабочего органа: 
1- рама; 2 -— консольная стойка; 3 рабочий орган 


Положение ненагруженного рабочего органа показано пунктиром. 
Оно соответствует проектной ориентации его в предположении, что упру- 
гие смещения отсутствуют. Подошва рабочего органа в этом случае должна 
лежать в горизонтальной плоскости. Ориентация рабочего органа задаётся 
вектором ориентации Л, расположенным в подошве рабочего органа. 

Под действием силы сопротивления Рь стойка деформируется и ра- 
бочий орган получает упругие смещения, которые можно представить как 
смещение $ носка лапы и поворота вектора ориентации на угол 6, из-за 
чего сила сопротивления Р изменяется по величине и направлению. Новое 
положение рабочего органа при небольших углах поворота удобно харак- 
теризовать тремя углами поворота относительно осей координат. Эти углы 
имеют традиционные названия: 6, - угол крена, ©, - угол тангажа, 6, - угол 
«рыскания». Они образуют вектор углов поворота 9 ={9х , бу, 02}' 

Смещения 1-й точки рабочего органа: 

5:=$+[9] 1, (1) 
где | - радиус-вектор текущей точки рабочего органа относительно его 
носка; [9] - матрица вращений на угол 6. 

Зависимость смещений от нагрузки $ (Р) характеризует траектории 
точек и в общем случае нелинейна. Для нелинейных систем справедливы 
законы деформирования в приращениях: 


4$ = [4] АР (2) 

аР = [С] 4, (3) 

40 = [Н] аР, (4) 

где [4] - симметричная матрица податливости размерности 3х3; 
[С]=[4]". - симметричная матрица жёсткости, размерности 3,3; 

[Н] - матрица угловой податливости размерности 3»3, причем В,=0, 


1=1,2,3. 
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Нелинейности делают все эти матрицы зависимыми от нагрузки. 
Для определённости их необходимо определять в рабочей точки вблизи 
номинальной нагрузки Р» ‚ при котором работает упругое крепление. 

Направление приложения силы описывает вектор направляющих 
косинусов п = {п:, п, пз} ", (п:= сои: , П2= СОЗ» , Пз= СОЗуз), Т.е. 

Р= [Р|п. (5) 

Упругие смещения характеризуются величиной и направлением. 
Направление упругих смещений на практике характеризуют параметром /. 
=$х/52 : 

^=(612С0$\/1+622 С05\у2+62з СО$\уз)/(611СО5\у1-+012С0у2+613С05\уз). = (6) 

В зависимости от направления смещения носка лапы существуют 
три кинематические режима: с заглублением ), > 0, равновесный ^, = 0 
и свыглублением /, < 0. В отношении устойчивости хода по глубине пред- 
почтительнее второй режим, не реагирующий на неоднородности почвен- 
ной структуры. 

Упругие смещения рабочего органа искажают геометрию резания и 
вызывают отклонения глубины хода от заданной (см. рис.1). Для обеспече- 
ния агротребований упругая кинематика должна регламентироваться. 

Агротехнические требования к кинематике упругих смещений рабо- 
чего органа сводятся к трём условиям: 

1. Глубина хода носка рабочего органа при рабочей нагрузке долж- 
на соответствовать заданной, т.е вертикальные смещения его должны рав- 
няться нулю: 457 = 0. 

2. Направление упругих смещений должно быть равно требуемому 
по технологии и энергосбережению, т.е. 

А, = Алреб . (7) 

3. Вертикальные упругие смещения любых точек подошвы УЗРО 
не должны выходить за пределы агродопуска Ад: 

{5 + [9] } ‚< Ад, (8) 
где © — матрица вращений для углов 0 = [Н]Р‚ ; Е - длина подошвы. 
Методика математического моделирования. В общем случае система 
«рабочий орган — почва — крепление» является системой с распределённы- 
ми параметрами и распределённой нагрузкой и должна описываться труд- 
норазрешимыми уравнениями в частных производных. Динамическую мо- 
дель теории терраупругости целесообразно строить переходом на макро- 
уровень через вектор состояния [4]. 

Вектор состояния системы «УРЗО-почва» включает векторы положе- 
ния носка рабочего органа П, вектор углов поворота 09 и их производные: 

$ =; 0; аП/аЕ; 46 /4%$}. (9) 

Сопротивление почвы описывает вектор Е размерности п=б, вклю- 
чающий составляющие силы сопротивления Р и момента М: 

Е = {Р ; М}={Рх, Ру, Р2, Мх, Му, М2}. (10) 

Согласно положениям теории разрушения почв в земледельческой 
механике [1, с.118-120], факторами, от которых зависит сопротивление 
почвы, являются глубина обработки, скорость и геометрия рабочего орга- 
на, площадь сечения пласта, свойства почвенной среды (Горячкин В.П., Ка- 
цыгин В.В., Зеленин А.Н., Новиков Ю.Ф., Бурченко П.Н. и др.). Режимные 
факторы входят в вектор состояния 3. Это означает, что для конкретной 
ситуации сопротивление почвы зависит от компонентов вектора состояния, 
являясь его вектор-функцией: 
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Е=ЕЗ]. (11) 
В существующей теории резания крепление считается абсолютно 
жёстким, обеспечивающим постоянные значения вектора состояния, т.е. $ 
о= {П; 0; №; 0} = соп$Е. Новым подходом к построению модели взаимодей- 
ствия рабочего органа с почвой в терраупругости является учёт изменений 
вектора состояния за счёт его упругих смещений, всегда возникающих 
под нагрузкой за счёт деформаций крепления. 
Упругие смещения $ и углы поворота 09 УЗРО образуют 6-мерный 
вектор обобщенных координат УЗРО: 
а ={5; @}= { х, у, 2, Ох, 6у@2}". (12) 
Вектор Ц изменения состояния рабочего органа за счет упругих 
смещений размерности 2п: 


У ={5; 0; 45/4; ае /аё}. (13) 
Результирующее значение вектора состояния 
$ =30 +1 = { Н+5; 0; М+ а$/аЕ; а9 /а4} (14) 


и сопротивление почвы движению рабочего органа зависит от упругих 
смещений: 
Е=Е {30 +1}. 

Использование вектора состояния позволяет при моделировании 
системы с распределёнными параметрами перейти на макроуровень. Разла- 
гая формы упругих смещений по формам собственных колебаний основных 
частот и приводя силы и массы к носку рабочего органа, удалось свести за- 
дачу к матричному дифференциальному уравнению в форме Коши [4], в 
котором за фазовые координаты приняты изменения вектора состояния Ц: 

УЦ =Сб9+8В{З} + р(®), (15) 
где с, \У- блочные матрицы собственных свойств упругой системы раз- 
мерности 2п со структурой вида 

Е 0 

и, 3 
-С -В 
где А - диагональная матрица коэффициентов инерции; В- симметричная 
матрица коэффициентов диссипации; С- симметричная матрица коэф- 
фициентов жесткости; Е, 0 - соответственно единичная и нулевая 
матрицы размерности п; В{3} - нелинейная детерминированная век- 
тор-функция реакции почвенной среды размерности 2п: 

В{3} = К{Зо + И}; (16) 
р(©) - случайная составляющая реакции вследствие неоднородностей 
почвенной среды. 

В структуре векторов реакции присутствуют векторы сопротивления 
почвы: детерминированная составляющая Е и случайная К ®: 

В{3} = [0; 2{3}], {Зо +4}=[0; {Зо +4}, р(®)=[0; К®]. 

Число п обобщенных координат выбирается в зависимости от 
конкретной задачи. При отсутствии моментов М упругие прогибы консоль- 
ного крепления пропорциональны углам поворота. Число независимых 
координат в этом случае п=3, и векторы обобщённых координат и сопро- 
тивления трёхмерны: 4=$=х, у, 2}'; Е=Р={Рх, Ру, Рх}". 

Система из уравнений (15) и (16) образует динамическую модель 
терраупругости. Первое уравнение системы описывает упругие смещения 
рабочего органа, второе - энергетику технологического процесса по прин- 
ципу "вход-выход". 


Е 0 
ол 


7! 





15 


Вестник ДГТУ, 2008. 1.8. №3(38) 








Структура модели представлена на рис 2. В её основе лежит 
встречно-параллельное соединение двух звеньев: звена УЗРО, описывае- 
мого уравнением (15), и звена энергетики технологического процесса поч- 
вообработки, описываемого уравнением (16). 

Как видно из схемы, упругость крепления образует обратную связь 
между выходом и входом модели. Наличие такой обратной связи является 
важной отличительной особенностью работы УЗРО, позволяющей объяс- 
нить все динамические эффекты процесса. 


Технологический процесс 
взаимодействия рабочего органа 
с ПОЧВОЙ 


Упругое крепление 
рабочего органа 





Рис 2. Структура модели терраупругости 


Анализ модели. При анализе нелинейную вектор-функцию В{3} под- 
вергли статистической линеаризации. По правилам статистической линеа- 
ризации 


В {3 } =В {30} + В, + [9] ч, (17) 
[9] = [Рш] + [Вы (18) 
а, = [9] 99, (19) 


где Ц - вектор математического ожидания фазовых координат; и - центри- 
рованный вектор фазовых координат; [0] - матрица коэффициентов 
статистической линеаризации порядка 2п 


0 0 
о ОО: 
[Ри] - матрица коэффициентов корреляции фазовых координат; [0..] 
- матрица взаимных коэффициентов корреляции упругой реакции и 
фазовых координат; [0], [0] - матрицы коэффициентов линеариза- 
ции порядка п; В, - вектор средних значений нелинейной вектор- 
функции В{Ц}; по физическому смыслу - реакция системы на упру- 
гость. 
Линеаризованные уравнения задачи включают динамическое уравнение 
колебаний “относительно среднего” и уравнение энергоэффекта “в сред- 


нем”: 
Уи" = {[6] + [9]}ч +р(©®, (20) 
[<] 9 +В {3+} + В, = 0. (21) 
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Достоинством модели линеаризованной модели терраупругости яв- 
ляется то, что каждое слагаемое имеет физический смысл: матрица [С] 
описывает свойства упругого крепления рабочего органа; матрица [0] 
описывает упругие и вязкие свойства почвенной среды; по физическому 
смыслу является матрицей квазиупругости почвенной среды; среднее зна- 
чение реакции системы на упругость В, характеризует энергоэффект УЗРО 
по отношению к жёсткому; В{Зо} - статическая характеристика жесткой 
системы, являющаяся традиционной задачей земледельческой механики о 
зависимости силы сопротивления от режимов резания. 

Влияние упругости крепления рабочего органа на энергозатраты 
процесса почвообработки, определяемые средним значением реакции на 
упругость В.возможно по выражениям (18), (19). Энергоэффект оказыва- 


ется зависимым не только от средних значений упругого сопротивления у, 


но и от их дисперсий В, что свидетельствует о существовании виброэф- 
фекта. 

Для практического использования уравнения модели (20), (21) мо- 
гут приводиться к трёхмерной форме, в которой параметры модели имеют 
традиционный физический смысл: 

[А]5” +[В|5' +[С5= рь(®, (22) 


[С] = [©] + [05] ‚ [Ву = [В] - [9 ; 
[С, [(] - матрицы жесткости нагруженной и ненагруженной системы; 
[Ва], [В] - матрицы затухания нагруженной и ненагруженной системы; 
[Ол, [03] - матрицы вносимого затухания и квазижесткости почвы; 
рв(©) - внешняя возмущающая сила с нулевым математическим ожи- 
данием ; 


где 


ре(Е) = рА(®)+ рш(®, 
где рл({) и рш(® - соответственно периодическая и случайная составляю- 
щая возмущений; $, $', 5” - центрированные значения соответственно 
вибросмещений, виброскорости и виброускорений. 
Возможности модели. Модель терраупругости, более полно описы- 
вая характер движения рабочего органа в почвенной среде, позволяет вы- 
явить ряд динамических эффектов [5]: прецессию собственных частот при 
нагружении системы, существование неустойчивых режимов движения ра- 
бочего органа типа дивергенции и флаттера, появление автоколебаний, их 
взаимодействие с вынужденными, и, что особенно ценно для практики, 
влияние упругости крепления рабочего органа на энергозатраты процесса 
почвообработки. 
Пример. Рассмотрим УЗРО культиватора в виде стрельчатой лапы. Сим- 
метрия нагружения позволяет ограничиться рассмотрением движения в 
продольно-вертикальной плоскости ХОУ с одним углом поворота бу - тан- 
гаж. Нагрузка на культиваторную лапу сводится к одной равнодействую- 
щей силе Р, крутящий момент отсутствует. Вектор упругих смещений будет 
состоять из двух координат $ = [х, 2]'. 
Всё это позволяет решать плоскую задачу, ограничиваясь двумя 
уравнениями в приращениях: 
а11 Х" + сих + С12 7 = ДАР/(х, 2)+ рьх(©, (23) 
а>2 7 “+ сах + С22 7 = ДР.(х, 7)-+ рьг(®). 
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По Горячкину В.П., тяговое сопротивление пропорционально пло- 
щади сечения деформируемого пласта 5п: 

Р = (К + =\?) бп , (24) 
где К, = - коэффициенты удельного сопротивления, зависящие от свойств 
почвы и геометрии резания. 

Обозначив (К + =\/2)=К(\), получим Р = К(\) 91. 

В свою очередь, площадь пласта и геометрия резания зависят от 
упругих смещений рабочего органа, тогда К = К(х, 7 ); $п =бп (х, 7) и сила 
сопротивления оказывается сложно зависимой от смещений х, 7: 

Р(х, 7) = К(СХ, 2) бп (х, 2). (25) 

Зададимся целью исследовать систему на устойчивость. Рассмотрим 
автономный режим, положив рьх({)= рь2(®)=0. При малых упругих смещени- 
ях (намного меньше глубины хода) нелинейностями вектор-функции Р(х, 7) 
можно пренебречь, тогда приращение составляющих силы сопротивления 
АР будет линейной комбинацией от х, 7: 

АР; = Ки Хх + К27, (26) 
ДР, = К2: Х + К 7. 

Подставляя эти выражения в уравнения (23) и пренебрегая силами 
демпфирования колебаний УЗРО в воздушной среде, получаем систему 
дифференциальных уравнений вида: 

Хх" + аих + ар 7= 0, (27) 
7" +алх + а 1 = 0. 
где 
11 = (С11- Ка1)/ ан: ; а 12 = (С12- К21)/ аз; 
21 = (С21- Ко1)/ а2з ; а 22 = (С>2- К22)/ аг2. 

Полученные уравнения - полный аналог классической модели 
флаттера в упрощённой форме по Пановко Я.Г. [7, с.334]. Такое совпаде- 
ние моделей носит принципиальный характер и свидетельствует о физиче- 
ском родстве процессов аэроупругости и терраупругости и обоснованности 
постановки задач терраупругости в почвообработке. 

Анализ корней характеристического уравнения системы (27) даёт 
два условия устойчивости [7, с.335]: 

ана» - ара»2: < (ан +а»2); /4, (28.1) 
а1: а? - а12 аз > 0. (28.2) 

Невыполнение первого условия, по Пановко Я.Г., приводит к дивер- 
генции, а невыполнение второго — к флаттеру [7, с.335]. 

Поскольку входящие в условия величины а к включают в себя зна- 
чения ск 1, К = 1, 2, то добиться соблюдения условий (28) устойчивой ра- 
боты УЗРО возможно только варьированием элементов матрицы жёсткости 
[С] , особенно её недиагональных элементов, характеризующих межкоор- 
динатные связи и кинематический режим работы крепления. 

Поиск периодических решений и энергоэффекта требует учёта не- 
линейностей функции Р(х, 7) по выражению (25). Аналитическое её опре- 
деление для практики имеет определённые трудности. 
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Целесообразнее получать Р(х, 7) экспериментально. При этом сле- 
дует учесть, что отношение 7/х = Л характеризует направление упругих 
смещений, и зависимость Р(х, 2) можно толковать как зависимость от кине- 
матического режима Р(х, Л). При горизонтальной силе сопротивления 
направляющие косинусы (5) п! =1, п›=Пз=0 и по выражению (6) показатель 
Л=0:12/ 911=- С12/С!1. Экспериментально изучена регрессионная модель вида: 

РЕЕ (Л, си1). (29) 

Получено регрессионное уравнение зависимости тягового сопротив- 

ления от кинематического режима и жёсткости: 
Р=68+8,0х!-7,16х!:?+2,5хх,-2,8х:2х› , (30) 
где хз=0,19%,; х.=(с!-45)/20. 

На рис.3 показаны зависимости Р(^.) для упругих подвесок разной 
жесткости: $ - си! = 25 Н/мм, ® - с:1 = 45 Н/мм и ДА -с:: = 65 Н/мм, при ско- 
рости движения \=12 км / ч, и глубине обработки Н =11 см. 

Важнейшим обстоятельством является наличие минимума силы со- 
противления вблизи нулевого кинематического режима. Минимум силы со- 
противления наблюдается при значениях кинематического показателя 
№= (-0;05)-(-0,1). Эти значения соответствуют слабому режиму с выглубле- 
нием. 
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Относительное тяговое сопротивление 


0,8 - 
ый ©` ©` ох ©` ©. < со хе 
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Рис.3. Зависимость тягового сопротивления 
от кинематического режима при различной жёсткости 
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Степень значимости оптимизации кинематического показателя /, ха- 
рактеризуется крутизной нарастания функции вблизи оптимума. Из графи- 
ков видно (рис.3), что изменение показателя кинематического режима от 0 
до 1,19 (усиление режима заглубления) вызывает значительное увеличе- 
ние тягового сопротивления от 10 до 45% в зависимости от жёсткости 
крепления. Если учесть, что для почвообработки существенным является 
снижение и на 5-10%, то отмеченные изменения говорят о высокой чув- 
ствительности тягового сопротивления к изменению кинематики упругих 
смещений и достаточно высокой эффективности способа энергосбережения 
за счёт её оптимизации. Чувствительным к кинематике упругих смещений 
оказывается и режим выглубления: изменение кинематического режима 
на 0,09 вызывает потери энергосбережения на 5 %. 

Таким образом, регрессионный анализ экспериментальных данных 

отчётливо подтверждает оптимизационность задачи: существование мини- 
мума тягового сопротивления при изменении упругой кинематики и высо- 
кую чувствительность энергетики процесса почвообработки к изменению 
кинематики упругих смещений, а значит, и потенциально высокую эффек- 
тивность способа энергосбережения за счёт её оптимизации. 
Выводы. Предлагаемая динамическая модель терраупругости отражает 
все особенности задачи взаимодействия рабочего органа с почвой (много- 
мерность, нелинейность, автоколебательность), пригодна для различных 
почвенных сред и различных рабочих органов; большинство существующих 
энергетических моделей почвообработки является её частными случаями. 

Применяемый для описания модели математический аппарат в виде 
матричных линеаризованных дифференциальных уравнений вводит её в 
русло современных достижений аналитической механики и, следовательно, 
может опираться на её развитое компьютерное программное обеспечение. 
Разработанная модель терраупругости в состоянии описать сложные дина- 
мические эффекты, сопровождающие движение УЗРО в почвенной среде: 
появление неустойчивых режимов, автоколебания, прецессию собственных 
частот и траекторий, вибро- и энергоэффект. 

Теория терраупругости позволила разработать практические реко- 
мендации, ряд которых положен в основу отраслевых методических указа- 
ний МУ 23.2.48-89, разработанных с НПО ВИСХОМ [6]. Они включают 
предложения, методики и алгоритмы по измерениям, расчётам, идентифи- 
кации и оптимизации упругих смещений рабочих органов почвообрабаты- 
вающих машин. 
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